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Abstract. 

We study models for surface growth with a wetting and a roughening transition 
using simple and pair mean- field approximations. The simple mean-field equations 
are solved exactly and they predict the roughening transition and the correct growth 
exponents in a region of the phase diagram. The pair mean-field equations, which are 
solved numerically, show a better accordance with numerical simulation and correctly 
predicts a growing interface with constant velocity at the moving phase. Also, when 
detailed balance is fulfilled, the pair mean field becomes the exact solution of the model. 
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1. Introduction 

Surface growing [TJ [2] is a phenomenon observed in nature as well as in the laboratory. 
As an example of the latter we cite the experimental technique known as molecular 
beam epitaxy which allows us the growing of surface at the atomic level. The surface 
is described by a height variable h(r, t) that gives the height of the deposited layer at a 
given point r of the substrate at time t. Several models have been introduced to describe 
the mean features of surface growing in which the heights are stochastic variables whose 
time evolution is governed by a Markovian stochastic process. Here we are concerned 
with the solution of such models by the use of mean-field approximations. 

Two relevant quantities are used to characterize the surface growth. One is the 
mean height h of the surface from the substrate and the other is the surface width w, 
which is a measure of the surface roughness. The divergence of w characterizes a rough 
interface. According to Family and Vicsek [3], the width of a rough surface of a sufficient 
large system behaves as 

u»~t 7 , (1) 

where 7 is the growth exponent. This behavior characterizes a rough thermodynamic 
phase. Otherwise, that is, if the width of an infinite system remains finite when t — > 00, 
the surface is smooth. A roughening transition takes place when, by varying the control 
parameters, the surface changes from a smooth to a rough surface. 

The interface may yet be pinned or moving. A moving thermodynamic phase is 
characterized by a constant velocity v of the interface, or in other words, by a linear 
growth of the mean height, that is, 

h = vt. (2) 

A depinning transition [4] from a pinned to a moving phase occurs when, by changing 
the control parameters, the interface begins to move with a constant velocity. This 
is also known as a wetting transition [5] since in the moving phase the mean height 
becomes infinitely large when t — > 00, that characterizes a wet phase. At the depinning 
transition, the mean height may not grow linearly with time but may behave according 
to 

h~t\ (3) 

In the model we study here the depinning transition coincide with the roughening 
transition. 

If we define P(h,t) as the one-point height probability distribution at time t, the 
mean height h and the square w 2 of the surface width are the first moment and the 
variance of this probability distribution, respectively. The asymptotic behavior given 
by Eqs. (JTJ and ([2]) can then be obtained from the following scaling form [6] 

p(M) = rv(^), (4) 

where f(x) is a universal function. The behaviors given by Eqs. ([I]) and ([3]) also follows 
from this same form by setting v = 0. 
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Some models have been introduced to describe the surface growth and the 
asymptotic behavior in which the heights are stochastic variables governed by Langevin 
equations. Two important ones are the Edwards- Wilkinson (EW) [7] 



where rj(r,t) is a white noise. Since the EW equation is linear, it stays invariant under 
the transformation h — > —h while the KPZ equation, which is nonlinear, lacks this 
property. In one dimension the growth exponent for the EW class is 7 = 1/4 and for 
the KPZ class is 7 = 1/3. 

In this work we study a growth model, introduced by Hinrichsen et al. |9j, with 
deposition and evaporation of particles that respects the restricted solid on solid (RSOS) 
condition [10] and in which the evaporation at the initial height is forbidden, that is, 
there is a wall at the zero height. The presence of the wall leads to a nonequilibrium 
wetting transition depending on the evaporation and deposition rates. It is worth 
mentioning that nonequilibrium wetting transition has been already studied by another 
approach namely by mapping the KPZ with a potential into a Langevin equation with 
a multiplicative noise [TH [12] . 

The growth model is studied here by means of simple and pair mean-field 
approximations. The mean-field approximation at the pair level, which does respect 
the RSOS condition exactly, is found to be capable of describing a moving phase, that 
is, a moving interface at constant velocity. Another feature of the two-site mean-field 
approximation is that it becomes the exact solution when detailed balance is fulfilled. 

Recently, a mean-field theory for surface growth has been introduced by Hinrichsen 
et al. [13] and used by Ginelli and Hinrichsen [6] to study a single step model for surface 
growth. The one-site and two-site mean-field approximations we use here are distinct 
from that of Hinrichsen et al. [TBj but share properties that are similar. 

The paper is organized as follows. In the next section we define the model to 
be examined and write down the master equation for general models with the RSOS 
condition. In section III we solve the master equation within the one-site mean-field 
approach exactly. The section IV is dedicated to the pair mean-field approach. The 
resultant equations are solved numerically and compared with the results obtained from 
the one-site approximation and from simulations. In section V we analyze the detailed 
balance condition and in section VI we present our conclusions. 

2. Model and master equation 

We consider a one-dimensional lattice where a discrete variable hi is attached to the site 
i, that represents the height of the pile of particles at site i. The stochastic variable hi 
takes the integer values. At time t = we consider a flat surface at the zero level, that 



_ = z/ V 2 /i + r / (r,t), 

and the Kardar-Parisi-Zhang (KPZ) [8] 
Bh 

= ^h + \{Vh) 2 + r,{v,t) 



(6) 



(5) 
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Figure 1. A configuration of heights with the RSOS condition 



is, hi = for all i. We consider models with random sequential updates of the heights, 
that have not only deposition but also evaporation. We will restrict ourselves only to 
models that obey the RSOS condition, that is, such that \hi — h i+ i\ < 1 (see Fig. [TJ). 
At each time step only one site is updated according to local stochastic rules. If site 
i is chosen, these rules will affect only the neighboring sites % + 1 and i — 1 and site i 
itself. The rate of the transition hi — > hi + n is denoted by c n (/tj_i, hi, h i+ i). The only 
possible transitions are those for which n — ±1 meaning that the height is increased 
(deposition) or decreased (evaporation) by just one unit. 

The possible transition rates are presented in Fig. [2] and are given by 

c+(h,h,h)=pi, (7) 

c + (h,h,h + l) = c+(h+l,h,h)=p3, (8) 

c+(h + l,h,h + l)=pB, (9) 

c.(h-l,h,h-l)=p2, (10) 

C-(h — 1, h, h) = C-(h, h, h — 1) = p 4 , (11) 

c-(h,h,h)=p e , (12) 

and the model so defined has six parameters: pi, ■ ■ ■ ,pe- By rescaling time we see that 
they are not all independent and one of them can set equal to unity. 

The probability P(h%, h%, . . . , t) of a given configuration (h%, h 2 , . . .) at time t obeys 
the master equation 

-P(/H,/l2,...,t) = 

at 

Yj J2i c n( h i~i, hi - n, h i+1 )P{h u . . . , h* - n, . . . , t) 

n=±l i 

- c n (/ij_i, hi, h i+ x)P{hx, ...,hi,...,t)}, (13) 

and we are considering periodic boundary conditions. For late use we write down the 
time evolution of the marginal probability distribution related to a given site, say site 
1. It is given by 

-n.P(hi, t) = i c n(h , hi - n, h 2 )P(h , hi - n, h 2 , t) 

dt h M n=±l 

- c n (h , hi, h 2 )P(h , hi, h 2 , t)}. (14) 
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Figure 2. Transition rates. 

We also write down the time evolution of the marginal probability distribution related 
to two consecutive sites, say site 1 and 2. It is given by 

i p {hi, h, t) = 2 {°n{ho, hi - n, h 2 )P(h , h x - n, h 2 , t) 

ho n=±l 

- Cn(h , h, h 2 )P(h , h x , h 2 , t)}, (15) 

valid for solutions that are translationally invariant, where we have made use of the 
symmetry c n (h , hi, h 2 ) = c n (h 2 , hi, h Q ). 

For certain values of the parameters the stochastic process exhibits detailed balance. 
This means to say that each term of the summation on the right hand side of Eq. ffT3l) 
vanishes in the time independent stationary state. In this case the time independent 
stationary probability distribution can be written as the product 

P{h x , h 2 , h 3 , h 4 , . . .) = i T(hi, h 2 ) T{h 2 , h 3 )T(h 3 , h A )... (16) 

where T{hi,h 2 ) are the elements of a symmetric matrix T, that vanish whenever the 
RSOS condition is not fulfilled. The detailed balance condition gives the following 
relations for the nonvanishing elements of T, denoted by = T{h, h) and = 
T(h,h + 1), 

p x Xl=p 2 Yl (17) 



p 3 X h Y h = p 4 Y h X h+ i, 

From these equations follows the relation [H] 

Pi Ps = pI_ 
P2 Pe pi ' 



(18) 
(19) 

(20) 
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which is the condition for detailed balance to hold. 

In the specific model [Qj we study here, the deposition occurs with a rate q and the 
evaporation, for nonzero heights, occurs with rates r or p, depending on the neighbors 
configuration. More specifically, if hi = 0, pi = p 3 = p 5 = q and p 2 = p$ = p& = 0; and 
if hi 0, 

Pi=P3=P5 = q, P2=Pi = r, p 6 = p. (21) 

The detailed balance condition (120]) gives p = r. Therefore, when p = r is respected the 
model can be solved exactly [9] for the pinned phase. 

Without the wall at the initial height, we would have a rough interface growing in 
the positive direction for q > q c and in the negative direction for q < q c . With the wall 
the phase q > q c is not affected, in the the sense that the interface still grows and is 
rough. For q < q c , in the presence of the wall, the interface is smooth and stays pinned 
to the substrate. Hence, at q = q c , the model displays a depinning and a roughening 
transition. The rough interface may display a EW or a KPZ behavior. In the case p = r, 
it is known that the crossover from EW to KPZ behavior [15] coincides with the critical 
point occurring at q — p. 

The case p = [161 EH] is special since in this case no particle can be evaporated 
from a completed filled layer and the presence of the wall thus makes no difference. 
The critical behavior at p = places the model in the universality class of the directed 
percolation (DP) class [TSlfTQ]. 

The order parameter of the pinned phase is the density of sites, Pq, in contact with 
the substrate, or the wall. It is assumed to behave near and below the transition as 

Po ~ (q c ~ qf- (22) 

The interface width w is finite at the pinned phase but diverges as one approaches the 
critical point. We assume that it diverges according to 

w ~ {q c - g)~ C . (23) 

The order parameter of the moving phase is chosen to be the velocity v of the interface 
growing, defined by 

. = >> = E*.5^.«). (24) 

which can also be written as 

v= nc n (ho,h 1 ,h 2 )P(h ,h 1 ,h 2 ,t). (25) 

hQ,hi,h,2 n=±l 

Near and above the transition we assume that the velocity behaves as 

v ~ (q - q c ) e . (26) 
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3. Simple mean-field approximation 

We are interested in studying the solution of the master equation associated to an infinite 
system. In order to solve the master equation (fT3]) we begin by using a simple mean 
field approximation. In this approach all the correlations are neglected resulting into the 
following approximation: P(hi-i,ht,hi+i,t) — P(hi-i,t)P(hi,t)P(hi+i,t). Insertion of 
this approximation into equation (fill) yields a closed equation for the one-site probability 
Pfc(i) distribution 

jP k = g(P|_ x - Pi) + (r - 2q)(P 2 P k+1 - Pl_ x P k ) 

+ (2r - q)(P k P 2 k+1 - P k -iP 2 k ) +p(P 3 k+1 - 7**?), (27) 

valid for the specific model defined by the rates (12TI) . where we are using the integer 
variable k in the place of the height h and initially Pq = 1 and P k = for k ^ because 
we are considering a initial flat surface. The variable 7fc = when k = and 7& = 1 
otherwise. The equation is valid for k > and for /c = provided we set P_i = on 
the right-hand side. 

Let us first look for a possible stationary state, that is, a time-independent solution, 
that correspond to a pinned phase. If we assume a solution of the form 

P k = AX\ (28) 

we can check easily by substitution that this is indeed a solution provided A is the root 
of the third-order algebraic equation 

-pX 3 + {q-2r)X 2 + {2q-r)X + q = 0. (29) 

The critical line, shown in Fig. [3j is found by letting A — > 1 with the result 

q c = \(p + Sr), (30) 



and the stationary solution occurs for q < q c . The normalization of P k , given by (J28J), 
gives A = 1 — A so that Po = 1 — A. The mean height h and the square of the interface 
width w 2 are determined as the average and the variance of the distribution given by 
(BED, that is 

h = (k) and w 2 = (k 2 ) - h 2 , (31) 

resulting in h — A/(l — A) and w = VX/(l — A). Near the critical line, A approaches 1 
as 1 — A = 2(q c — q)/{r + p) from which follows the results P ~ (q c — q), h ~ (q c — q)^ 1 
and w ~ (q c — q)' 1 - 

We solve the equation for q > q c by a method similar to that of Ginelli and 
Hinrichsen [H] in which a continuous version of the equation is used. Writing h = k, the 
equation for the probability distribution P(h,t) becomes, up to second order, 

2 



op n2 dP 




p ^ + -p 



1 „od 2 P 



dh 2 



(32) 
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Figure 3. Phase diagram for r = 1 in the p — q plane. The critical lines, separating 
the moving (M) and pinned (P) phases, were determined by simple mean- field (SMF), 
pair mean- field (PMF) approximation and numerical simulations (S). 



where 



a = 12(q - q c ), 



(33) 



and 



b = 2q + r + 3p. (34) 
At the critical line, q = q c , the coefficient a vanishes and we end up with the 



equation 



dP 
~dt 



P 



'dP" 
dh 



1 d*P 

2 dh 2 



(35) 



which can be solved by assuming the scaling form (J3J) for P(h,t). A consistency is 
achieved only if we choose the growth exponent as being 7 = 1/4. Along the critical 
line, the simple mean-field approximation is then compatible with an EW behavior. The 
substitution of the scaling form into the equation leads to the following equation for the 
scaling function f(x) 

f(x) + xf{x) + Abf(x)[f'(x)] 2 + 2b[f(x)] 2 f"(x) = 0, (36) 

whose solution is 



/(*) 



1- 



x" 
2b' 



(37) 



The asymptotic probability distribution, valid for large times, can then be written as 

(38 



P(h,t) 



I 1 h 2 
t^'Zbt 
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For q > q c it suffices to consider only the linear term in b. The equation for P(h, t) 
reads 

^ - -aP^ (39) 

dt ~ ^ dh' [ } 

where now a is strictly positive. This equation can be solved by assuming again a scaling 

relation of the form (j3J). Now, however, the consistency requires an exponent 7 = 1/3 

so that the simple mean-field approximation is then compatible with a KPZ behavior. 

The substitution gives the following equation for the scaling function 

f(x) + xf'(x) = 3a[f(x)] 2 f'(x), (40) 

whose solution is 

/(*) = vf ■ ( 41 ) 

The asymptotic probability distribution, valid for large time, can then be written 

as 

P(h,t) = M. (42) 
V at 

We remark that P(h, t) vanishes for values of h larger than a certain h max which depends 
on time. This maximum value of h is determined by the normalization of P(h,t) 
and is given by h max = (9at/4) 1//3 . This allows us to determine h = 3/i max /5 and 
w = (12/175) 1/2 /i max so that 

h~w~ (q-q c ) 1/3 t 1/3 . (43) 

Within the simple mean-field approximation we were able to observe a roughening 
transition occurring at the transition line shown in Fig. [3J Nevertheless, this 
approximation is not able to predict a moving surface with a constant velocity. This 
deficiency will be overcome in the next section by means of the pair mean-field 
approximation. 

4. Pair mean-field approximation 

In the pair mean-field approach we solve the equation f|T5|) by using the approximation 
for the three-site probability distribution 

P MM = ^^M , (44) 

where P(h\, h2,t) and P(h,2,t) are the two-site and one-site probability distribution. 
These two quantities are related by 

P(h 2 , t) = J2 *) = £ P(h 2 , h, t). (45) 

Due to the RSOS condition there are actually three types of two-site probabilities: 
Pk,k{t), Pk,k-i(t) and Pk,k+i{t), where as before we are using the integer variable k in 
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Figure 4. Order parameter Pq as a function of q for p = 0.5 obtained from simple 
mean- field (SMF), pair mean-field (PMF) approximation and numerical simulations 
(S) performed with 2048 sites and 100 independent realizations. The critical values 
are q c = 0.875 for SMF, q c = 0.8539 for PMF and q c = 0.8346(1) for the simulations. 



the place of the height h. Together with the one-site probability Pk they form a set of 
four variables. However they are not all independent since Pk+i,k — Pk,k+i and 

Pk = Pk,k + Pk,k-l + Pk,k+l- (46) 

Therefore we are left with two independent variable, for each k, which we choose to 



be Pk t k and Pk,k+i- Initially we have P; 



'k,k+i = f° r all k and P^k = for k ^ and 
P ,o = I- 

To simplify notation, we denote Pk,k+i by yk and P^k by Xk so that 

Pk = Xk + Vk + Vk-i- (47) 

The pair mean-field equations for Xk and yk, corresponding to the model defined by 
(EH), then reads 



d 



x k = 2q 



dt 
+ 2r 



y k -i{vk-i + Xk-i) x k (x k + y k ) 



P 



k-l 



yk{yk + x k +i) x k yk-i 



Pi 



where we should set x_ 



k+l 



y-i 



P k 

when k 



0, and 



d 



:Vk = g- 



xi 



Vk 



— r 



Vk 



+ p- 



x 



k+l 



(4? 



(49) 



dt Pk Pfc+i Pfc+i 

The numerical integration of the coupled equations (j48|) and (j49|) . performed by repeated 
iteration from the initial condition x k = for k ^ 0, xq = 1 and yk = for all fc, allows us 
to determine the one-site probability distribution (|47|) . From the obtained distribution 
we get the mean height and the interface width by the use of Eqs. (|3T|) and the interface 
velocity by derivating numerically the mean height with respect to time. 
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Figure 5. Interface width w as a function of q forp — 0.5 obtained from pair mean-field 
(PMF) approximation and numerical simulations (S) performed with 2048 sites and 
100 independent realization. The critical values are q c — 0.8539 and q c — 0.8346(1), 
respectively. 



The transition line is shown in Fig. [31 For comparison we also show the transition 
line obtained from numerical simulation. Below the transition line, in the pinned phase, 
the density of sites P$ in contact with the substrate is found to behave as 

P ~(g c -g), (50) 

so that the exponent f3 = 1, the same value found in the one-site mean-field 
approximation. The width of the interface diverges at the critical line as 

w~(q c -q)- l '\ (51) 

yielding an exponent £ = 1/3. 

As an example, we show in Figs. H] and [5j respectively, the order parameter P and 
the width w as functions of the deposition rate q, at p = 0.5. We see that P is linear 
near the transition point and that w~ 3 is linear supporting the behavior ( l5Tl) . 

Above and at the transition line, the numerical integration gives a time dependent 
probability distribution Pk{t) from which we calculate the average height and the width 
of the interface. Along the transition line we found that 

w ~ t 1/4 , q = q„ (52) 

and above it 

w ~ t 1/3 , q > q c . (53) 

In Figs. [6] and [7] we show a data collapse of Pfc(i) by using the scaling form (j3j) at the 
critical point (p = 1 and q = 1) and above the critical point, inside the moving phase 
(p = 1 and q = 2). In the former case we have used 7 = 1/4 and in the latter, 7 = 1/3. 
The good data collapse leads us to the results (1521) and ( 1531) . 
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0.12 



0.08 - 



0.04 




Figure 6. Scaling function / = P(h, t) t 1 / 4 versus x — ht^ 1 / 4 as obtained from the 
pair mean-field approximation for the values of t shown in figure, at the critical point 
q = p = 1. The function f(x) extrapolates to zero when x — > 0. 

0.12 | . 1 . 1 . 1 



0.08 



0.04 - 



-15 




Figure 7. Scaling function / = t) t 1 / 3 versus x — (h — vt)t~ 4 ^ as obtained 
from the pair mean-field approximation for the values of t shown in figure, inside the 
moving phase at the point q = 2 and p = 1. The velocity is found to be v = 0.112. 



The asymptotic results for the width of the interface are the same as those found 
by means of the one-site mean-field approximation. However, as already mentioned, 
the two-site approximation is capable of describing a moving interface. We found that, 
inside the moving phase the mean height increases with a constant velocity v so that 
the pair mean-field approximations predicts correctly a linear growth of the moving 
interface. As we approach the critical line the velocity vanishes according to 

u~(?-9c), (54) 

that is, 0=1. 
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As we have seen, the pair mean-field approximation takes into account the 
correlation of two consecutive sites which results in the exact satisfaction of the RSOS 
condition in opposition to the simple mean-field approximation. Because of that the 
pair mean-field results are closer to the simulation ones when compared to the simple 
mean-field results. For instance, the transition line is in general closer to the simulation 
results and has the right concavity. 



5. Detailed balance 



When detailed balance is obeyed, it turns out that the two-site approximation becomes 
an exact solution. In this case the stationary probability distribution is of the form 
ffT6]) and the elements of X k and of the transfer matrix T are related to Xk and 
by Xk = XfcPfc/A and yk = Yk\J PkPk+\l A, where A is the dominant eigenvalue of T. 
The detailed balance conditions (IT71) . (jTSJ) and fTTTJT) are then equivalent to the following 
relations: 

p = r = l, q— = - — , y k =x k x k+ i 1 (55) 

from which follows the solution 

— = -q k , Vk = -q h+1/2 , (56) 

Pk A y r PkPk+i A 

where A is a constant to be found. That this is indeed a solution of the pair mean-field 
equations can be checked by substitution. 

Inserting these relations into equation ( 1471) we get the following eigenvalue equation 
for (p k = s/T k 

q k - 1/2 <t>k-i + q k <Pk + q k+1/2 <Pk+i = A <p k (57) 

so that the constant A is identified as the eigenvalue, actually, the dominant eigenvalue 
A. The dominant eigenfunction gives the one-site probability distribution Pk = 0|. The 
eigenvalue equation ( 1571) was found by Hinrichsen et al. [13] who solved it in the vicinity 
of the critical point by using a continuous height approximation. Their solution gives 
the following results in the vicinity of the critical point q = q c = 1 

w ~ (q c ~ <?r V3 , (58) 

and 

P{0)~(q c -q)- (59) 

These results are the same results found in the previous section by numerical integration 
of the two-site pair approximation equations. 



6. Conclusion 



We have studied a lattice model for surface growth by means of one-site and two- 
site mean-field approximations. Within the simple mean-field approximation, whose 
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equations were solved exactly, we were able to observe a roughening transition occurring 
with growth exponents equal to 1/4 at the transition line and 1/3 inside the rough 
phase. This simple approximation was not able to predict a moving surface with a 
constant velocity. However, this deficiency was overcome by means of the pair mean- 
field approximation which predicts correctly a linear growth of the moving interface. 
In general, the pair approximation gives results that are in better accordance with 
numerical simulations. The critical line is closer to that given by numerical simulations. 
The growth exponent is 1/4 at the critical line and 1/3 inside the moving phase. At 
p = 1 this approximations correctly predicts the crossover from the EW, at q = 1, to 
KPZ, for q > 1. 

The pair mean-field approximations gives also the following results. As one 
approaches the critical line from inside the pinned phase, the order parameter vanishes 
with an exponent (3 = 1 and the width of the interface diverges with an exponent 
£ = 1/3. If one approaches the line from the moving phase, the velocity of the interface 
vanishes with an exponent — 1. 

The pair mean-field approximation has two important features: it takes into account 
the RSOS condition exactly and reduces to the exact solution when detailed balance is 
fulfilled. This last feature indicates that the pair mean-field is a good approximation 
for the nonequilibrium case (p ^ 1), where exact results are generally not known. This 
approximation can also be use to get the properties and the phase diagram of the general 
model defined in Fig. [21 with or without a wall. For growth models that do not respect 
the RSOS condition, in principle, the pair approximation can also be used but the 
equations will have a more cumbersome form. 
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